1 


1 


1 


\ 


NASA TECHNICAL 
MEMORANDUM 


JASA TM X- 73329 


MEASUREMENTS OF THE STOCHASTIC NATURE OF 
ATMOSPHERIC SPECTRAL AMPLITUDES 


By William C. Cliff, George H. Flchtl, 

Margaret Alexander, and Salvador Arias 

Space Sciences Laboratory • 

\ 

i 


August 1976 


NASA 



George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 

(NA5A-TK- X-73329) MEASUREMENTS OF THE N76-31783 

STOCHASTIC NATURE OF ATMOSPHERIC SPECTRAL 
AMPLITUDES (NASA) 33 p HC $4.00 CSCL 04A 


MS K • Form 1|40 iK'r \ f i me 10 71 ) 


Unclda 

G 3/46 02431 


1 


TECHNICAL REPORT STANDARD TITLE PAGE 


1. REPORT HO. 

NASA TM X-73329 


4 TITLE AND SUBTITLE 


2. GOVERNMENT ACCESSION NO. 


Measurements of the Stochastic Nature of Atmospheric 
Spectral Amplitudes 


7 . author ( s) William C. Cliff, George H. Fichtl, Margaret 




9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


12 SPONSORING AGENCY NAME AND ADDRESS 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


15. SUPPLEMENTARY NOTES 


3. RECIPIENT'S CATALOG NO. 


5. REPORT DATE 

August 1976 


6. PERFORMING ORGANIZATION CODE 


8. PERFORMING ORGANIZATION REPORT # 


10. WORK UNIT NO. 

1 1. CONTRACT OR GRANT NO. 

13. TYPE OF REPOR'i 8c PERIOD COVERED 

Technical Memorandum 


14. SPONSORING AGENCY CODE 


Prepared by Space Sciences Laboratory, Science and Engineering 


16. ABSTRACT 

The stochastic nature of the power spectral amplitudes of the neutral atmospheric 
boundary layer is examined. Probability density distributions and probability distributions 
of longitudinal and lateral power spectral amplitudes are computed from neutral atmospheric 
boundary layers at Cape Kennedy. The statistical distributions are computed for frequencies 
of 0,006, 0.01, 0.03, 0.06, 0.1, and 0.5 Hz at each of the elevations of 18, 30, 60, 90, 

120, and 150 m. When the probability density distributions are properly nondimensionalized, 
the data tend to collapse to a universal curve. An empirical curve fit to the universal 
nondimensionalized probability density distribution is also given. Probability distributions 
of individual frequency power spectral amplitudes are also presented for all elevations and 
frequencies. 

An interesting observation from the data is that greater than 10 percent of the time 
the power spectral amplitude at a given frequency will generally be more than three times 
the temporal mean value computed by standard Fourier techniques. The standard pow r er 
spectral density curves are also included in the report. 



18. DISTRIBUTION STATEMENT 


C 

Unclassified — Unlimited 


19. SECURITY CLASSIF. (of thU report! 

20 SECURITY CL AS 51 F, (of thta page) 

21. NO. OF PAGE 7 S 

22 PRICE 

Unclassified 

Unclassified 

33 

NTIS 


MSFC - Form 3 212 (Rev December 1173 ) 

> ' T- 


For sale by National Technical Information Service, Springfield, Virginia 11 15! 

1 I I 




















X 


TABLE OF CONTENTS 


Page 


INTRODUCTION 1 

PROCEDURE 2 

RESULTS 3 

CONCLUSIONS 7 

REFERENCES 25 


iii 



LIST OF ILLUSTRATIONS 


i 

sl 

I 

$ 

Figure 

1 . 

2 . 

i 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 


Title 

Power spectral density of longitudinal wind component 
(data taken at an elevation of 18 m, rms = 1. 89 m/s) 

Probability density distribution of the longitudinal power 
spectral amplitudes for a frequency of 0. 06 Hz (data 
taken at an elevation of 1 8 m) 

Probability distribution of the longitudinal power spectral 
amplitudes for a frequency of 0.06 Hz (data taken at an 
elevation of 18 m) 

Comparison of probability density distribution of longitudinal 
power spectral amplitudes for various frequencies (data 
taken at an elevation of 18 m) 

Comparison of 1. 0 -probability density distribution versus 
amplitude for various frequencies (data are of the longi- 
tudinal wind component taken at an elevation of 18 m) 

Comparison of the probability density distribution of the raw 
data amplitudes for various frequencies (data are of the 
longitudinal wind component taken at an elevation of 18 m). . 

Probability density distribution of longitudinal power 
spectral amplitudes for various frequencies (data taken 
at an elevation of 30 m) 

Probability density distribution of longitudinal power 
spectral amplitudes for various frequencies ( data taken 
at an elevation of 60 m) 

Probability density distribution of longitudinal power 
spectral amplitudes for various frequencies ( data taken 
at an elevation of 90 m) 

Probability density distribution of longitudinal power 
spectral amplitudes for various frequencies (data taken 
at an elevation of 120 m) 




Page 

9 

9 

10 

10 

11 

11 

12 

12 

13 

13 


iv 





LIST OF ILLUSTRATIONS (Continued) 


Figure Title Page 

11. Probability density distribution of longitudinal power 
spectral amplitudes for various frequencies ( data taken 

at an elevation of 150 m) 14 

12. Probability density distribution of lateral power spectral 
amplitudes for various frequencies ( data taken at an 

elevation of 18 m) 14 

13. Probability density distribution of lateral power spectral 
amplitudes for various frequencies ( data taken at an 

elevation of 30 m) 15 

14. Probability density distribution of lateral power spectral 
amplitudes for various frequencies ( data taken at an 

elevation of 60 m) 15 

15. Probability density distribution of lateral power spectral 
amplitudes for various frequencies (data taken at an 

elevation of 90 m) 16 

16. Probability density distribution of lateral power spectral 
amplitudes for various frequencies ( data taken at an 

elevation of 120 m) 16 

17. Probability density distribution of lateral power spectral 
amplitudes for various frequencies (data taken at an 

elevation of 150 m) 17 

18. Power spectral density of longitudinal wind component 

(data taken at an elevation of 30 m, rms = 1. 93 m/s) ] 7 

19. Power spectral density of longitudinal component (data 

taken at an elevation of 60 m, rms = 1. 79 m/s) 18 

20. Power spectral density of longitudinal component (data 

taken at an elevation of 90 m, rms = 1. 90 m/s) 18 

21. Power spectral density of longitudinal component (data 

taken at an elevation of 120 m, rms = 1.66 m/s) 19 


v 



LIST OF ILLUSTRATIONS (Concluded) 


st 


Figure 

Title 

Page 

22. 

Power spectral density of longitudinal component (data 
taken at an elevation of 150 m, rms= 1.63 m/s) 

. 19 

23. 

Power spectral density of lateral component (data taken 
at an elevation of 18 m, rms =1.48 m/s) 

. 20 

24. 

Power spectral density of lateral component (data taken 
at an elevation of 30 m, rms - 1.48 m/s) 

20 

25. 

Power spectral density of lateral component (data taken 
at an elevation of 60 m, rms - 1. 35 m/s) 

21 

26. 

Power spectral density of lateral component (data taken 
at an elevation of 90 m, rms = 1.48 m/s) 

21 

27. 

Power spectral density of lateral component (data taken 
at an elevation of 120 m, rms =- 1. 38 m/s) 

22 

28. 

Power spectral density of lateral component (data taken 
at an elevation of 150 m, rms = 1.46 m/s) 

22 

29. 

Probability density distribution of power spectral 
amplitudes 

23 

30. 

Distribution of variance as a function of amplitude ...... 

23 

31. 

Distribution of mean spectral value as a function of 
amplitude 

24 

32. 

Empirical probability distribution of power spectral 
amplitudes . 

24 


i 

\ 


f 


vi 



TECHNICAL MEMORANDUM X-73329 




MEASUREMENTS OF THE STOCHASTIC NATURE OF 
ATMOSPHERIC SPECTRAL AMPLITUDES 

INTRODUCTION 


In 1915, G. I. Taylor [1] first examined eddy motion in the atmosphere, 
and in 1917 [2] he performed experiments to map the shape of eddies in the 
atmosphere. In 1938, Taylor [3] published the relationship between correlation 
measurements and spectra (harmonic analyses) to statistically analyze an 
atmospheric velocity record. Since that time, harmonic analysis has been the 
standard technique to statistically compress velocity records. The relationship 
of the harmonic amplitudes of the velocity record to that of the harmonic ampli- 
tudes representing the autocorrelation is given in other works [4] . The objec- 
tive of this report is to specifically address the stochastic nature of these har- 
monic amplitudes (i. e. , the stochastic nature of the power spectral amplitudes 
of wind records) . The general statistical technique for analyzing a wind record 
is to take a wind velocity record and perform the standard Fourier decomposition 
of the data record given by: 


R(t) = t / u ( t )u(t+r)dt (1) 

0 

and 

R (t) = 7 ^ / C(w) e 1WT 6w 

uu' ' 2n J ' ’ 

-GO 

where 


= f R ( t) e ^ T dr 
J uu 

-CO 


C(o>) 


( 2 ) 
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and where R (t) is the time averaged autocorrelation of the u component of 
uu' ' 

velocity and C(o> ) is the average power spectra! amplitude at a frequency, w . 

It is well understood, however, that even under neutral atmospheric conditions 
the value of statistical measures of the atmospheric winds are stochastic and 
thus depend to some degree upon the integration time or observation time from 
which the statistic is evaluated. Thus, by allowing the power spectral amplitude 
to be a function of time as well as frequency, one is able to develop the statistical 
properties of C(co ,t) such that 


c(w) = h f c(w,t) dt . (3) 

0 


PROCEDURE 


The technique used here to evaluate C(u> ,t) was to first select long time 
histories of neutral atmospheric wind velocities. These time histories were 
filtered at a desired frequency; the filtered function was then squared and divided 
by twice the filter bandwidth. The peaks of the resulting time function now 
correspond to a temporal record of the power spectral amplitudes of the original 
data record; that is, the peaks represent C(u.' ,t), where 


C(cj,t) 


A 2 (c o ,0 
2AB 


and where A (a> ,t) is the temporal Fourier amplitude of the wind record, such 
that 


c<°) = 7 / C(co.t) dt = i / dt . (4) 

0 0 

As expected, equation (4) simply states that the general power spectral 
amplitude is equal to the time averaged stochastic power spectral amplitude. 
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The power spectral amplitude, C(w ,t), may be represented by the well known 
relationship between power spectral amplitude and its frequency complement in 
Fourier representation of the original velocity time history. 


RESULTS 

Using the technique outlined in the previous section, the probability 
density distribution (PDD) of the power spectra for individual frequency com- 
ponents are computed from temporal wind records of five neutral atmospheric 
boundary layers. Longitudinal and lateral wind components are analyzed at 
frequencies of 0.006, 0.01, 0.03, 0.06, 0. 1, and 0. 5 Hz from data records 
taken at elevations of 18, 30, 60, 90, 120, and 150 m. All corresponding 
standard statistical measurements are also included. A typical example of the 
calculated longitudinal wind component power spectra is presented in Figure 1. 

The main objective of this report, as noted previously, is to examine 
the stochastic nature of the power spectra amplitude for individual frequency 
components. The stochastic nature of the individual frequency spectral ampli- 
tudes is presented in the form of PDD. The PDD of an individual frequency is 
presented in Figure 2. Figure 2 is the PDD of the longitudinal wind component 
for a frequency of 0.06 Hz and the data taken at an elevation of 18 m. Figure 2 
is taken from the same data record as Figure 1. Thus, the mean value 
of the power spectral amplitude from Figure 2, 


00 

C(co ) = C(w,t) = f C(w,t) P(C(w,t)) dC(w,t) , 

-00 


must be equivalent, to the ensemble mean value (Fig. 1) which was calculated by 
conventional harmonic analysis. The two values are seen to be equivalent. It 
is interesting to note that because the power spectra amplitude is always positive 
(which results in a natural skewing) , the rms value of the amplitudes can easily 
be larger than the mean value. For the case in Figure 2, the mean value is 
approximately 10, whereas the rms value of the same curve has a value greater 
than 13. Because the PDD of the power spectra components is always positive, 
it is tempting to speculate that the PDD could be approximated by standard 
statistical distributions such as the gamma, exponential, or beta. An empirical 
curve fit to the data taken is presented later. The density distribution appears 
to have a negative exponential form. 
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Figure 3 presents the probability distribution ( PD) for the same data as 
used to compute Figures 1 and 2. An interesting observation from Figure 3 is 
that greater than 10 percent of the probability curve lies beyond a value three 
times that of the mean value. 

Similarity techniques are used here to collapse the PDD of all frequencies 
(of a given wind direction and elevation) to a single curve. The collapse of the 
normalized PDD for the power spectral amplitude of the longitudinal component 
for several frequencies taken at an elevation of 18 m is presented in Figure 4, 
while an example of the collapse of 1. 0 minus the PD for the same data is given 
in Figure 5. 

The PDD and PD of the raw data amplitudes for individual frequencies 
and the PDD and PD of the power spectral amplitudes for individual frequencies 
were computed. These amplitudes appeared to be symmetrical with the PDD 
and PD of the negative amplitudes which were assumed to be mirror images of 
the positive PDD and PD; thus, only positive amplitudes were computed. An 
example of the PDD of the raw data amplitudes for various frequencies is 
presented in Figure 6. 

Plots of the normalized PDD of the frequency-dependent longitudinal 
power spectral amplitudes of data taken at elevations 30, GO, 90, 120, and 150 m 
are presented in Figures 7 through 11, respectively. Plots of the- normalized 
PDD of the frequency-dependent lateral power spectra, amplitudes of data taken 
at elevations of 18, 30, 60, 90, 120, and 150 m are presented in Figures 12 
through 17, respectively. 

Figures 18 through 22 present the standard power spectral density curves 

of the longitudinal wind component for elevations 30, 60, 90, 120, and 150 m, 

respectively. Note that the value taken from these curves at a given frequency 

is the mean value of the stochastic power spectral amplitude and is denoted 

C(w) . 
m 


Figures 23 through 28 present the standard power spectral density curves 

of the lateral wind component for elevations of 18, 30, 60, 90, 120, and 150 m, 

respectively. Again, as was previously the case, a value taken from these curves 

represents the mean value of the stochastic power spectral amplitude which is 

denoted as C(u>) . 

' 'm 

The most interesting observation from this collection of data is 
that all of the normalized probability density curves of the power spectral ampli- 
tudes (whether of the longitudinal or lateral component, and independent of the 
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elevation at which the data were taken and independent of the frequency of con* 
cem) appear to collapse to a universal curve which may be expressed as 


PlC(u.t)) . 


P(C(u.) m ) 


e- 11 "' 8 * 1 + 30.25 e* 20 * 


= 0 


for x > 0 
for x < 0 


( 5 ) 


where 


x = C{g> t ,t) 

C(w) 

' 'm 


and where P(C(o> ,t)) is the value of the PDD for a given amplitude C(w ,t) , 
p(C(o) ) ) is the value of the PDD of the average power spectral density, 

C(w ,t) is the given amplitude, and C(o> ) is the average amplitude taken 

directly from the standard power spectral density curve. 

Equation (5) then leads to formulation of an empirical PDD for 
C(u> ,t) /C(u> ) by simply dividii j equation (5) by the area under the curv r e. 
Thus, m 


, /c(q> ,t) \ _ 


P(x) = 


_x°.8+l -20x 

e + 30,25 e 

4.592 


= 0 


for 

X ^0 (6) 

for 

X 

A 

O 


where 


X = 


Cloull 

C<w) 


m 


5 


It should be noted that 


) 

I 

j 


CO 


/ P ( x ) dx = l 


-00 


and also that 


X = / x P(x) dx = l 


-CO 


The second equation results from the requirement that 


00 

7 / C(u,t) = C(u) 

-CO 


Equation (6) conforms to both of these requirements. A plot of equation 
(6) is presented in Figure 29, and tabular results arc given in the Table. Other 
curves which may interest a designer are curves representing the distribution 
of the variance of the spectral amplitudes { i. e. , X 2 P(x) J and the distribution 
of the mean spectral values [ i. e. , x P(x) | . These curves are presented in 
Figures 30 and 31, respectively, to show what values of spectral amplitude con- 
tribute to the value and mean value of spectral amplitudes. Another obvious 
curve of interest to the designer is the probability distribution curve, that is, 
the integral of the probability density distribution. The empirical probability 
distribution curve resulting from equation (6) is presented graphically in 
Figure 32 and in tabular form in the Table. From the Table it is interesting to 
note that approximately 9 percent of the fluctuations are greater than three times 
the mean spectral amplitude. This value is found by looking at X = 0 and the 

x 

corresponding value of j P(x) dx from the Table. 
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CONCLUSIONS 
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The stochastic nature of power spectral amplitudes for neutral atmos- 
pheric boundary layers has been studied. " 2 study shows that the PDD of the 
spectral amplitudes monotonically decreases for increasing amplitude. Because 
of the natural skewing of the PDD, the rms value (standard deviation) may be 
larger than the mean value. Approximately 10 percent of the expected instan- 
taneous power spectral amplitudes has been shown to exceed a value greater 
than three times the mean value. Perhaps the most significant contribution is 
that the PDD of the spectral amplitudes will collapse to a universal curve. This 
curve is independent of whether the longitudinal or lateral wind component is 
analyzed, independent of the specific frequency analyzed, and independent of the 
elevation at which the atmospheric wind record was taken. 
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TABLE. VALUES OF THE PROBABILITY AND PROBABILITY DENSITY 
DISTRIBUTION OF PO.VER SPECTRAL AMPLITUDE,* 



C (cj,t) 

X- — ; 

c 

X- JxP Wdx*1 ; 
0 

P(x)-e x °- 8+1 +30.25e- 20 * 

4.592 

JP (x) d x * 1 
o 




X 




X 

X 

PM 

/ P(x)dx 


X 

PM 

/ P(x)dx 



O 




o 

.10 

1.41096 

.33902 


3.60 

.03686 

.93670 

.20 

.57564 

.42516 


3.70 

.03466 

.94024 

.30 

.42473 

.47297 


3 J 0 

.03259 

.94356 

.40 

.37207 

.51213 


3.90 

.03066 

.94669 

.50 

.33699 

.54716 


4.00 

.02885 

.94964 

JO 

.30769 

.57903 


4.10 

.02716 

.95241 

.70 

.28196 

.60819 


4.20 

.02557 

.95502 

.90 

.25904 

.63495 


4.30 

.02406 

.95747 

.90 

.23849 

.65955 


4.40 

.02269 

.95979 

1.00 

.21997 

.68223 


450 

.02138 

.96197 

1.10 

.20322 

.70316 


4.60 

.02015 

.96402 

1.20 

.18801 

.72251 


4.70 

.01900 

.96596 

1.30 

.17416 

.74043 


4.80 

.01792 

.96779 

1.40 

.16152 

.75703 


4.90 

.01691 

.96951 

1.50 

.14996 

.77244 


5.00 

.01595 

.97114 

1.90 

.13936 

.78675 


5.10 

.01505 

.97267 

1.70 

.12963 

.80006 


5.20 

.01421 

.97412 

1 J 0 

.12068 

.81244 


5.30 

.01342 

.97549 

1.90 

.11244 

.82398 


5.40 

.01267 

.97678 

2.00 

.10484 

.83472 


5.50 

.01197 

.97800 

2.10 

.09782 

.84475 


5.60 

.01131 

.97915 

2.20 

.09133 

.85411 


5.70 

.01069 

.98023 

2.30 

.08533 

.86285 


5.80 

.01010 

.98126 

2.40 

.07976 

.87102 


5.90 

.00955 

.96224 

2 J 0 

.07460 

.87865 


6.00 

.00903 

.98315 

2.90 

.06981 

.88580 


6.10 

.00854 

.98402 

2.70 

.06537 

.89249 


6.20 

.00807 

.98486 

2 J 0 

.06123 

.89875 


6.30 

.00764 

.98562 

2.90 

.05739 

.90462 


6.40 

.00723 

.98636 

3.00 

.06380 

.91012 


6.50 

.00684 

.98706 

3.10 

.06047 

.91528 


6.60 

.00647 

.98771 

3.20 

.04736 

.92012 


6.70 

.00613 

.98834 

3.30 

.04446 

.92466 


6.80 

.00580 

.98893 

3.40 

.04175 

.92893 


6.90 

.00550 

.96949 

3.50 

.03922 

.93293 


7.00 

.00521 

.99002 
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POWER SPECTRAL DENSITY 
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Figure 1. Power spectral density of longitudinal wind component 
( data taken at an elevation of 18 m, rms = 1. 89 m/ s) . 



AMPUTUOC 

Figure 2. Probability density distribution of the longitudinal power 
spectral amplitudes for a frequency of 0. 06 Hz 
(data taken at an elevation of 18 m). 
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Figure 3. Probability distribution of the longitudinal power spectral 
amplitudes for a frequency of 0.06 Hz (data taken at an 
elevation of 18 m) . 
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Figure 4. Comparison of probability density distribution of longitudinal 
power spectral amplitudes for various frequencies 
(data taken at an elevation of 18 m). 
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Figure 5. Comparison of 1. 0 probability density distribution versus 
amplitude for various frequencies (data are of the longitudinal wind 
component taken at an elevation of 18 m) . 



NORMALIZED SPECTRAL DENSITY (A/AB1 / IA/AB) M 


Figure 6. Comparison of the probability density distribution of the raw 
data amplitudes for various frequencies (data are of the longitudinal 
wind component taken at an elevation of 18 m) . 
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Figure 7. Probability density distribution of longitudinal power spectral 
amplitudes for various frequencies ( data taken at an elevation of 30 m) . 
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Figure 8. Probability density distribution of longitudinal power spectral 
amplitudes for various frequencies (data taken at an elevation of 60 m). 






NORMAL POWER SPECTRAL DENSITY AMPLITUDE [C (u.tJ/C M J 

Figure 9. Probability density distribution of longitudinal power spectral 
amplitudes for various frequencies ( data taken at an elevation of 90 m) . 



Figure 10, Probability density distribution of longitudinal power spectral 
amplitudes for various frequencies (data taken at an elevation of 120 m) . 
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Figure 11. Probability density distribution of longitudinal power spectral 
amplitudes for various frequencies (data taken at an elevation of 150 m) . 
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Figure 12. Probability density distribution of lateral power spectral amplitudes 
for various frequencies (data taken at an elevation of 18 m) . 
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Figure 13. Probability density distribution of lateral power spectral amplitudes 
for various frequencies (data taken at an elevation of 30 m) . 



Figure 14. Probability density distribution of lateral power spectral amplitudes 
for various frequencies (data taken at an elevation of 60 m) . 
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Figure 15. Probability density distribution of lateral power spectral amplitudes 
for various frequencies (data taken at an elevation of 90 m) . 
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Figure 16. Probability density distribution of lateral power spectral amplitudes 
for various frequencies (data taken at an elevation of 120 m) . 
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Figure 17. Probability density distribution of lateral power spectral 
amplitudes for various frequencies (data taken at 
an elevation of 150 m). 



Figure 18. Power spectral density of longitudinal wind component 
(data taken at an elevation of 30 m, rms = 1. 93 m/s) . 
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Figure 19. Power spectral density of longitudinal component (data 
taken at an elevation of 60 m, rms = 1. 79 m/s) . 
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Figure 20. Power spectral density of longitudinal component (data 
taken at an elevation of 90 m, rms = 1.90 m/s) . 
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Figure 21, Power spectral density of longitudinal component (data 
taken at an elevation of 120 m, rms = 1, 66 m/s) . 



Figure 22, Power spectral density of longitudinal component (data 
taken at an elevation of 150 m, rms = 1.63 m/s). 
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Figure 23. Power spectral density of lateral component (data 
taken at an elevation of 18 m, rms = 1. 48 m/ s) . 



Figure 24. Power Bpectral density of lateral component (data 
taken at an elevation of 30 m, rms = 1.48 m/s). 
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Figure 25. Power spectral density of lateral component (data 
taken at an elevation of 60 m, rms «* 1. 35 m/s) . 



Figure 26. Power spectral density of lateral component (data 
taken at an elevation of 90 m, rms = 1.48 m/s). 
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! 27. Power spectral density of lateral component (data 
taken at an elevation of 120 m, rms = 1. 38 m/s) . 
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Figure 28. Power spectral density of lateral component ( data 
taken at an elevation of 150 m, rms = 1.46 m/s) . 
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